Macroscopic law of conservation revealed in the population dynamics of Toll-like receptor signaling

Stimulating the receptors of a single cell generates stochastic intracellular signaling. The fluctuating response has been attributed to the low abundance of signaling molecules and the spatio-temporal effects of diffusion and crowding. At population level, however, cells are able to execute well-defined deterministic biological processes such as growth, division, differentiation and immune response. These data reflect biology as a system possessing microscopic and macroscopic dynamics. This commentary discusses the average population response of the Toll-like receptor (TLR) 3 and 4 signaling. Without requiring detailed experimental data, linear response equations together with the fundamental law of information conservation have been used to decipher novel network features such as unknown intermediates, processes and cross-talk mechanisms. For single cell response, however, such simplicity seems far from reality. Thus, as observed in any other complex systems, biology can be considered to possess order and disorder, inheriting a mixture of predictable population level and unpredictable single cell outcomes

Emergence of macroscopic simplicity from the Tumor Necrosis Factor-alpha signaling dynamics

The Tumor Necrosis Factor-α (TNF-α), a cytokine produced during the innate immune response to invading pathogens, is involved in numerous fundamental cellular processes. Here, to understand the temporal activation profiles of the TNF-α regulated signaling network, we developed a dynamic computational model based on the perturbation-response approach and the law of information (signaling flux) conservation. Our simulations show that the temporal average population response of the TNF-α stimulated transcription factors NF-κB and AP-1, and 3 groups of 180 downstream gene expressions follow first-order equations. Using the model, in contrast to a well-known previous study, our model suggests that the continuous activation of the third group of genes is not mainly due to the poor rate of mRNA decay process, rather, the law of signaling flux conservation stipulates the presence of secondary signaling, such as feedback mechanism or autocrine signaling, is crucial. Although the living system is perceived as sophisticated and complex, notably, our work reveals the presence of simple governing principles in cell population dynamics

Biophysical mechanism for Ras-nanocluster formation and signaling in plasma membrane

Ras GTPases are lipid-anchored G proteins which play a fundamental role in cell signaling processes. Electron micrographs of immunogold-labeled Ras have shown that membrane-bound Ras molecules segregate into nanocluster domains. Several models have been developed in attempts to obtain quantitative descriptions of nanocluster formation, but all have relied on assumptions such as a constant, expression-level independent ratio of Ras in clusters to Ras monomers (cluster/monomer ratio). However, this assumption is inconsistent with the law of mass action. Here, we present a biophysical model of Ras clustering based on short-range attraction and long-range repulsion between Ras molecules in the membrane. To test this model, we performed Monte Carlo simulations and compared statistical clustering properties with experimental data. We find that we can recover the experimentally-observed clustering across a range of Ras expression levels, without assuming a constant cluster/monomer ratio or the existence of lipid rafts. In addition, our model makes predictions about the signaling properties of Ras nanoclusters in support of the idea that Ras nanoclusters act as an analog-digital-analog converter for high fidelity signaling.Comment: 8 figures. PLoS ONE, in pres

Genetic vehicle comprising majority of lowly expressed genes guides cell fate decision

Cells remarkably take a specific differentiation path among the multiple possibilities that can arise through the multi-dimensional regulation of genome activities. Such deterministic processes suggest the existence of cellular attractors. However, the origins and drivers of the attractors still remain elusive. Here we analyzed the temporal neutrophil differentiation microarray data for two different stimuli, dimethyl sulfoxide (DMSO) and all-trans-retinoic acid (atRA), and expressed their collective dynamics by temporal Pearson correlation and mutual information coordinates. We constructed ensemble of the genes which showed reduction of correlation fluctuations following the inverse square root law. Evaluating their temporal probability density distributions resulted in the emergence of distinct high density localizations from non-localized low density spread distributions, forming attractor cores for both atRA and DMSO. These attractor cores overlapped, pointing to the existence of a neutrophil cell fate attractor. Notably, we found the localizations of correlation distributions for the majority of lowly expressed genes (LEGs) ensembles overlapped with the whole genome attractor cores, while the remaining genomes' localizations did not overlap. Therefore, we postulate the existence of genetic vehicle, made up mainly of LEGs, for the guidance of cell fate

ABioTrans: A Biostatistical Tool for Transcriptomics Analysis

Here we report a bio-statistical/informatics tool, ABioTrans, developed in R for gene expression analysis. The tool allows the user to directly read RNA-Seq data files deposited in the Gene Expression Omnibus or GEO database. Operated using any web browser application, ABioTrans provides easy options for multiple statistical distribution fitting, Pearson and Spearman rank correlations, PCA, k-means and hierarchical clustering, differential expression (DE) analysis, Shannon entropy and noise (square of coefficient of variation) analyses, as well as Gene ontology classifications

Analysis of the Molecular Networks in Androgen Dependent and Independent Prostate Cancer Revealed Fragile and Robust Subsystems

Androgen ablation therapy is currently the primary treatment for metastatic prostate cancer. Unfortunately, in nearly all cases, androgen ablation fails to permanently arrest cancer progression. As androgens like testosterone are withdrawn, prostate cancer cells lose their androgen sensitivity and begin to proliferate without hormone growth factors. In this study, we constructed and analyzed a mathematical model of the integration between hormone growth factor signaling, androgen receptor activation, and the expression of cyclin D and Prostate-Specific Antigen in human LNCaP prostate adenocarcinoma cells. The objective of the study was to investigate which signaling systems were important in the loss of androgen dependence. The model was formulated as a set of ordinary differential equations which described 212 species and 384 interactions, including both the mRNA and protein levels for key species. An ensemble approach was chosen to constrain model parameters and to estimate the impact of parametric uncertainty on model predictions. Model parameters were identified using 14 steady-state and dynamic LNCaP data sets taken from literature sources. Alterations in the rate of Prostatic Acid Phosphatase expression was sufficient to capture varying levels of androgen dependence. Analysis of the model provided insight into the importance of network components as a function of androgen dependence. The importance of androgen receptor availability and the MAPK/Akt signaling axes was independent of androgen status. Interestingly, androgen receptor availability was important even in androgen-independent LNCaP cells. Translation became progressively more important in androgen-independent LNCaP cells. Further analysis suggested a positive synergy between the MAPK and Akt signaling axes and the translation of key proliferative markers like cyclin D in androgen-independent cells. Taken together, the results support the targeting of both the Akt and MAPK pathways. Moreover, the analysis suggested that direct targeting of the translational machinery, specifically eIF4E, could be efficacious in androgen-independent prostate cancers

Integrated Transcriptome and Binding Sites Analysis Implicates E2F in the Regulation of Self-Renewal in Human Pluripotent Stem Cells

Rapid cellular growth and multiplication, limited replicative senescence, calibrated sensitivity to apoptosis, and a capacity to differentiate into almost any cell type are major properties that underline the self-renewal capabilities of human pluripotent stem cells (hPSCs). We developed an integrated bioinformatics pipeline to understand the gene regulation and functions involved in maintaining such self-renewal properties of hPSCs compared to matched fibroblasts. An initial genome-wide screening of transcription factor activity using in silico binding-site and gene expression microarray data newly identified E2F as one of major candidate factors, revealing their significant regulation of the transcriptome. This is underscored by an elevated level of its transcription factor activity and expression in all tested pluripotent stem cell lines. Subsequent analysis of functional gene groups demonstrated the importance of the TFs to self-renewal in the pluripotency-coupled context; E2F directly targets the global signaling (e.g. self-renewal associated WNT and FGF pathways) and metabolic network (e.g. energy generation pathways, molecular transports and fatty acid metabolism) to promote its canonical functions that are driving the self-renewal of hPSCs. In addition, we proposed a core self-renewal module of regulatory interplay between E2F and, WNT and FGF pathways in these cells. Thus, we conclude that E2F plays a significant role in influencing the self-renewal capabilities of hPSCs

Emergent Genome-Wide Control in Wildtype and Genetically Mutated Lipopolysaccarides-Stimulated Macrophages

Large-scale gene expression studies have mainly focused on highly expressed and ‘discriminatory’ genes to decipher key regulatory processes. Biological responses are consequence of the concerted action of gene regulatory network, thus, limiting our attention to genes having the most significant variations is insufficient for a thorough understanding of emergent whole genome response. Here we comprehensively analyzed the temporal oligonucleotide microarray data of lipopolysaccharide (LPS) stimulated macrophages in 4 genotypes; wildtype, Myeloid Differentiation factor 88 (MyD88) knockout (KO), TIR-domain-containing adapter-inducing interferon-β (TRIF) KO and MyD88/TRIF double KO (DKO). Pearson correlations computed on the whole genome expression between different genotypes are extremely high (>0.98), indicating a strong co-regulation of the entire expression network. Further correlation analyses reveal genome-wide response is biphasic, i) acute-stochastic mode consisting of small number of sharply induced immune-related genes and ii) collective mode consisting of majority of weakly induced genes of diverse cellular processes which collectively adjust their expression level. Notably, temporal correlations of a small number of randomly selected genes from collective mode show scalability. Furthermore, in collective mode, the transition from large scatter in expression distributions for single ORFs to smooth linear lines emerges as an organizing principle when grouping of 50 ORFs and above. With this emergent behavior, the role of MyD88, TRIF and novel MyD88, TRIF-independent processes for gene induction can be linearly superposed to decipher quantitative whole genome differential control of transcriptional and mRNA decay machineries. Our work demonstrates genome-wide co-regulated responses subsequent to specific innate immune stimulus which have been largely neglected

Analysing Dynamical Behavior of Cellular Networks via Stochastic Bifurcations

The dynamical structure of genetic networks determines the occurrence of various biological mechanisms, such as cellular differentiation. However, the question of how cellular diversity evolves in relation to the inherent stochasticity and intercellular communication remains still to be understood. Here, we define a concept of stochastic bifurcations suitable to investigate the dynamical structure of genetic networks, and show that under stochastic influence, the expression of given proteins of interest is defined via the probability distribution of the phase variable, representing one of the genes constituting the system. Moreover, we show that under changing stochastic conditions, the probabilities of expressing certain concentration values are different, leading to different functionality of the cells, and thus to differentiation of the cells in the various types